Low-temperature magnetotransport of narrow-gap semiconductor FeSb 2 
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We present a study of the magnetoresistance and Hall effect in the narrow-gap semiconductor FeSb2 at low 
temperatures. Both the electrical and Hall resistivities show unusual magnetic field dependence in the low- 
temperature range where a large Seebeck coefficient was observed. By applying a two-carrier model, we find 
that the carrier concentration decreases from 1 down to 10~ 4 ppm/unit cell and the mobility increases from 2000 
to 28000 cm 2 /Vs with decreasing temperature from 30 down to 4 K. At lower temperatures, the magnetoresis- 
tive behavior drastically changes and a negative magnetoresistance is observed at 3 K. These low-temperature 
behaviors are reminiscent of the low-temperature magnetotransport observed in doped semiconductors such as 
As-doped Ge, which is well described by a weak-localization picture. We argue a detailed electronic structure 
in FeSb2 inferred from our observations. 

PACS numbers: 72.20.My, 72.20.Pa 



I. INTRODUCTION 

Fe-based narrow-gap semiconductors have attracted con- 
siderable attention because of their unusual transport and 
magnetic behaviors, which closely resemble those of sev- 
eral rare-earth compounds termed "Kondo insulators". In 
4 /-electron Kondo insulators such as YbB 12 , Ce 3 Bi 4 Pt 3 , 
and CeNiSnji— a narrow gap is formed by the hybridiza- 
tion between localized 4/ and conduction electrons, leading 
to unique physical properties. The magnetic susceptibility 
obeys the Curie- Weiss law described by the local 4/ moment 
at high temperatures, followed by a broad maximum, and is 
suppressed at low temperatures. A large Seebeck coefficient 
due to a large density of states near the Fermi level was ob- 
served in the low-temperature range. In analogy to the 4f- 
electron system, the narrow-gap semiconductor FeSi has been 
intensively studied as a prototype of the 3<i-electron Kondo 
insulator. 1 This material displays an unusual crossover from 
high-temperature metal to low-temperature insulator, 6 ^ rem- 
iniscent of the 4/-electron Kondo insulators. In the photoe- 
mission experiments, however, the measured electronic struc- 
ture has no distinct features relevant to a Kondo picture and 
is qualitatively explained within the band calculations by the 
density functional theory without many-body effects, 8 £ puz- 
zling the precise role of electron correlation in the formation 
of a narrow gap in this compound. 

FeSb2 is another candidate for the Fe-based Kondo 
insulator, 10 the crystal structure of which is shown in the inset 
of Fig. 1. This material shows a unique temperature depen- 
dence of the magnetic susceptibility similar to that of FeSi. 
The resistivity is well characterized by thermal activation with 
two energy gaps of Ai ^30 meV and A2 ~ 5 meVj^— The 
novelty of this compound is best highlighted by the colossal 
Seebeck coefficient S ~ —45 mV/K near 10 K, where the 
smaller gap is open. 1 1 Such a huge value of S immediately 
indicates a novel mechanism, such as the strong electron cor- 
relation expected in Kondo insulators, which is pointed out 
by the thermodynamic and optical measurements ! 12 ' 15,16 The 
Seebeck coefficient and resistivity, however, are extremely 
sensitive to the sample quality, and the reported maximum 
S values span from -500 pN/K up to -45 mV/K^ri^i 7 . 



A systematic study of the impurity effects on the transport 
properties of FeSb2 revealed that the Seebeck coefficients are 
simply related to the carrier concentration doped by impuri- 
ties, not to the electron correlations, as seen in conventional 
semiconductors. 18 

In this paper, we present the magnetotransport study in 
FeSb2 single crystals, especially for the investigation of its un- 
usual electronic state at low temperatures. We conduct a two- 
carrier analysis for the conductivity tensor, which provides us 
with a proper evaluation of the mobility and carrier concentra- 
tion. The carrier concentration decreases from 1 down to 10 -4 
ppm/unit cell with decreasing temperature from 30 down to 4 
K. The small gap of 6 meV obtained from the temperature 
variation of the carrier concentration agrees well with the pre- 
vious results. The mobility reaches a large value of 28000 
cm 2 /Vs at 4 K, which enables us to observe the large See- 
beck coefficients even in such a low carrier concentration in 
FeSb2. At lower temperatures, the magnetoresistive behavior 
dramatically changes from positive to negative. The negative 
magnetoresistance is generally found in doped semiconduc- 
tors at low temperatures and derives from a weak localiza- 
tion of carriers, 19 ' 20 implying that a similar electronic struc- 
ture with a considerable impurity level is responsible for the 
unusual transport properties in FeSb2. 

II. EXPERIMENT 

High-quality single crystals of FeSb2 were grown by a self- 
flux method using metal powders of 99.999 % (5N) pure Fe 
and 99.9999 % (6N) pure Sb, as described in Ref.Ql The 
single-crystalline nature was checked by Bragg spots in the 
Laue pattern. The resistivity was measured from 100 down 
to 2 K. The magnetic field dependence of electrical and Hall 
resistivities were measured up to 70 kOe in the tempera- 
ture range from 30 down to 3 K. These transport properties 
were measured using a conventional four-probe dc method in 
the Quantum Design Physical Property Measurement System 
(Quantum Design, Inc.). For the magnetoresistance measure- 
ments, we employed the transverse configuration, i.e., with 
the current perpendicular to the magnetic field. We observed 
no significant direction dependence with respect to the crys- 
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FIG. 1 : (Color online) The temperature dependence of the resistivity 
p (circle symbols, left axis) and the Seebeck coefficient S (triangle 
symbols, right axis) under zero magnetic field. In the inset, the crys- 
tal structure of FeSb2 is schematically shown. 
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FIG. 2: (Color online) (a), (b) The transverse magnetoresistance 
Ap/p, and (c), (d) the Hall resistivity p xy , as a function of mag- 
netic field measured at various temperatures. The vertical axes of (a) 
and (c) are nearly two orders of magnitude larger than those of (b) 
and (d), respectively. 



talline axes. For measurements of the Hall resistivity, the mag- 
netic field was swept from +70 to —70 kOe, and the Hall volt- 
age was determined as Vh = \ [V (+H) — V{—H)\, where V is 
the voltage across the Hall terminals. The Seebeck coefficient 
was measured using a steady-state method. 



III. RESULTS AND DISCUSSION 

Figure 1 shows the temperature dependence of the resistiv- 
ity p and the Seebeck coefficient S below 100 K under zero 
magnetic field. The resistivity exhibits an insulating behavior, 
including two upturns separated by a plateau near 20 K, indi- 
cating the existence of two energy gaps in FeSb2, as shown 
in previous reports*^— The Seebeck coefficient changes its 
sign around 30 K and a large value of l^l ~ 1400 /iV/K is 
observed near 20 K. 

In Figs. 2(a) and 2(b), we show the transverse magnetore- 
sistance Ap/p = [p xx (H ) - p xx (0)] jp xx (0), where p xx is the 
dc resistivity, measured at various temperatures from 3 to 30 
K. Above 15 K, a conventional quadratic magnetic field de- 
pendence of Ap/ p is observed in low fields, but in high fields, 
it increases almost linearly with increasing magnetic field. At 
8 and 10 K, Ap/p tends to saturate in high fields and reaches 
1.5 around 70 kOe. Such large values of Ap/ p indicate high 
carrier mobility in this system. At 4 and 5 K, Ap/p is strongly 
suppressed at high magnetic fields. At 3 K, a negative magne- 
toresistance is observed. 

The magnetic field dependence of the Hall resistivity p xy 
at several temperatures is shown in Figs. 2(c) and 2(d). The 
Hall resistivity is highly nonlinear in H from 4 to 20 K at 
the moderately high magnetic fields, in contrast to the typical 
i/-linear dependence of p xy confirmed at 30 K. 

To properly evaluate the carrier concentration n and the 
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FIG. 3: (Color online) The magnetic field dependence of the conduc- 
tivity tensor (a) a xx and (b) a xy from 8 to 30 K. Solid curves are the 
calculation using Eqs. (T} and ([Ill- 



mobility p, we now discuss the magnetic field dependence 
of the conductivity tensor. In Figs. 3(a) and 3(b), we show 
Oxx = Pxx/{p 2 xx + ply) and a xy = -p xy /{p 2 xx + p 2 xy ), 
respectively, in the temperature range between 8 and 30 K. 
The magnetic field dependence of a xx and a xy systematically 
changes with decreasing temperature. Figure 4 shows a xx and 
a xy below 5 K. The magnetic field dependence of a xx below 
4 K is clearly different from that of the above temperatures, 
and a drastic change of a xy curves is observed from 4 to 3 K. 
In the Boltzmann transport theory, the magnetic field depen- 
dence of the conductivity tensor is expressed by using n and 
p. We here analyze a xx and a xy by using a two-carrier model 
in which one carrier is of high mobility and the other is of low 
mobility (pH <C 1). In this case, the H dependence of the 
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FIG. 4: (Color online) The magnetic field dependence of the conduc- 
tivity tensor (a)-(c) a xx and (d)-(f) a xy from 3 to 5 K. Solid curves 
are the fitting calculation using Eqs. l[T) and l(2j. a xx below 4 K and 
a xy at 3 K cannot be fitted with these equations. 



conductivity tensor is described as 21 



a xy (H) = n xy en 2 H 



1 + (p xx H) 2 



1 + {^xyHf 



C x 



Cxy 



(i) 



(2) 



Here 

^ xx foxy), P xx ([t X y\ and C xx {C X y} ^re the carrier con- 
centrations, the carrier mobilities, and the low-mobility com- 
ponents for a xx (<J xy ), respectively. The best agreements of 
Oxx and a xy with Eqs. (Q]) and (O, respectively, are shown by 
solid curves in Figs. 3 and 4. However, the data of <j xx at 3 
and 4 K and those of a xy at 3 K are unadaptable to the above 
two-carrier model, and we will discuss this reason later. 

In Figs. 5(a) and 5(b), we show the l/T dependence of 
the carrier mobility and the carrier concentration obtained 
from the fitting for the conductivity tensor, respectively. The 
carrier mobilities and fi xy , which are obtained indepen- 
dently from G xx and a xy , are in excellent agreement with each 
other, justifying our evaluation for these quantities within the 
two-carrier model. This is also evidenced by the quantitative 
agreement between n xx and n xy . The low-mobility compo- 
nents C xx and C xy from the fitting results are displayed in 
Fig. 5(c). The mobility /ii ow estimated from 

Cx X and Cxy is 

about 25 % of the high-mobility components at 10 K. Thus, 
Hiov/H *C 1 holds at low fields and even at high fields, and 
our fitting model using Eqs. (1) and (2) is still valid because 
the contributions of the low-mobility components to a xx and 
a xy are much smaller than those of high-mobility carriers. 
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FIG. 5: (Color online) l/T dependence of (a) the carrier mobility fi 
and (b) the carrier concentration n, and (c) the low-mobility compo- 
nents C' xx and C X y, obtained from the fitting results for the conduc- 
tivity tensor. The dotted line (b) is the fitting result by the thermal 
activation function n oc exp(— A/2/cbT), which yields A = 6 meV. 



Let us first discuss the carrier mobility in this material. 
As shown in Fig. 5(a), the mobility of the dominant car- 
riers is about 2000 cm 2 /Vs at 30 K and increases with de- 
creasing temperature. At 4 K, the mobility reaches a sig- 
nificantly large value of p ~ 28000 cm 2 /V s, comparable to 
that of high-purity Si or Ge semiconductors* 2 ^ 2 ^ Such a large 
value of the mobility at low temperatures has been also ob- 
served in several Kondo insulators and charge-density-wave 
systems, 24,25 because of an enhancement of the scattering time 
due to gap opening. Previous studies of the transport proper- 
ties have shown a large mobility in FeSb2»^& The Hall mo- 
bility fin = \Rh\/ p, which is obtained from the Hall coeffi- 
cient Rh and the resistivity p,— is different from the present 
results: hh exhibits a broad peak around 8 K and its max- 
imum value is five times smaller than ours. This deviation 
from our results is clearly due to the evaluation method in 
which only a single band is assumed. Hu et al. reported a gi- 
ant carrier mobility of [i ~ 8 x 10 4 cm 2 /Vs near 10 K, which 
was estimated only from the magnetic field dependence of the 
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Hall coefficient within a two-carrier model.- Their raw data 
of p xy (H) at 10 K is, however, much smaller than ours at the 
same temperature and can be plotted into our data between 30 
and 25 K, which implies that \i is of the order of 10 3 cm 2 /Vs. 
The fitting model used in Ref. 26 to evaluate both the carrier 
concentration and the mobility is only valid for very-high- 
mobility electron systems in a nonquantizing low-field range, 
in which p xy (H) exhibits a strong nonlinear H dependence, 27 
incompatible with the nearly iJ-linear dependence of p xy (H) 
observed in Ref. 26. Moreover, in the high-field range, the 
model equation of the Hall coefficient does not depend on H, 
and the values of the mobility and carrier concentration can- 
not be extracted independently using this model. 27 Thus, the 
model might greatly overestimate the carrier mobility in this 
system. 

We next discuss the carrier concentration shown in Fig. 
5(b). With decreasing temperature from 30 down to 4 K, 
it rapidly decreases from 10 16 cm~ 3 (1 ppm/unit cell) to 
10 12 cm -3 (10~ 4 ppm/unit cell). As shown by the dot- 
ted line, it is well fitted by the thermal activation function 
n oc exp(— A/2fcsT), which yields A = 6 meV. The sys- 
tematic study of the impurity effects on the transport proper- 
ties indicates that the electric conduction in this temperature 
range is governed by extrinsic carriers; the main carriers are 
excited from the impurity level, which locates slightly below 
the bottom of the conduction band. 18 In this picture, A corre- 
sponds to the energy gap between the impurity level and the 
bottom of the conduction band. In other words, the gap A is 
unlikely to be a gap in the Kondo insulators. 

Below 4 K, we could not fit a xx using Eq. ([T). As seen in 
Fig. 1(b), the magnitude of Ap/p is extremely small com- 
pared to those above 5 K, and the negative magnetoresis- 
tance is observed around low magnetic fields at 3 K. These 
behaviors have been widely observed in doped semiconduc- 
tors in the low-temperature range because of a weak elec- 
tronic localizations^^ . In the case of As-doped Ge, impu- 
rity carriers show an electric conduction at very low temper- 
atures, where the thermal excitation into the conduction band 
is negligibly small. The carriers are weakly localized by a 
random potential, leading to a negative magnetoresistance. 20 
The striking similarity of low-temperature magnetoresistive 
behaviors between FeSb2 and conventional semiconductors 
suggests a considerable impurity level in the electronic struc- 
ture of FeSb2: The impurity level exists slightly below (~ 6 
meV) the bottom of the conduction band and carriers in the 
impurity level are predominant on the electrical transport at 3 
K. 

Finally, we discuss the Seebeck coefficient S with a calcu- 
lation using a nondegenerate model expressed as 28 

where 6 is a scattering parameter. Now we calculate the value 
of S using A = 6 meV, which is estimated from the \IT de- 
pendence of the carrier concentration, and obtain \S\ ~ 600 
/iV/K at 20 K. Here we assume 5 = 3/2, which is used in 
conventional semiconductors. The calculated value is close 
to, but nearly two times smaller than, the experimental results 



of S = -800 1400 pV/K in Refs. M and 0. An 

electron-electron correlation is unlikely to be the origin of this 
deviation because it has an insignificant effect on the Seebeck 
coefficient in an insulating regime. 30 On the other hand, the 
phonon-drag effect may play an essential role for an enhance- 
ment of the Seebeck coefficient in the low-temperature re- 
gion, as seen in many pure semiconductors^^ The phonon- 
drag term in the Seebeck coefficient can be approximately 
expressed as S p = 0v p l p /[iT, where v p and l p are the ve- 
locity and the mean free path of a phonon, respectively, and 
< j3 < 1 is a parameter which characterizes the relative 
strength of the electron-phonon interaction! 3 ^! From the high 
thermal conductivity below 30 K, a long mean free path of 
l p ~ 10 /im at 20 K was estimated by Bentien et a h 1 21 In 
addition, they calculated the sound velocity to be 3000 m/s. 
Using these parameters, S p at 20 K is estimated to be nearly 
3 mV/K at maximum, which suggests that the phonon-drag 
effect can yield a measurable contribution to the Seebeck co- 
efficient in this material. On the other hand, it is still diffi- 
cult to understand the colossal Seebeck coefficient {S ~ —45 
mV/K at 10 K) that was observed in Ref. 1 1 within the above 
description. 



IV. CONCLUSION 

We have measured the magnetic field dependence of the 
electrical and Hall resistivity of pure FeSb2 single crystals and 
analyzed the conductivity tensor by using a two-carrier model 
to evaluate the carrier concentration and mobility. The car- 
rier concentration decreases from 1 down to 10~ 4 ppm/unit 
cell with decreasing temperature from 30 down to 4 K. The 
mobility significantly increases with decreasing temperature 
and reaches 28000 cm 2 /Vs at 4 K. At lower temperatures, the 
transport behavior changes drastically and a negative magne- 
toresistance is observed. In doped semiconductors, the nega- 
tive magnetoresistance due to a weak localization is generally 
found at low temperatures, implying a considerable impurity 
level in the electronic structure of FeSb2. The temperature 
variation of the carrier concentration indicates the existence 
of an energy gap (~ 6 meV) between the bottom of the con- 
duction band and the impurity level. The magnetotransport re- 
ported in this paper strongly suggests that the low-temperature 
transport of FeSb2 is well understood as an extrinsic semicon- 
ductor with ppm-level impurity. The phonon-drag effect may 
work, but the electron correlation effect is secondary. 



Acknowledgements 

The authors would like to thank M. Sato for the initial mo- 
tivation and collaboration at an early stage of this work. The 
authors also appreciate H. Kontani and A. Kobayashi for their 
helpful discussion. They are also indebted to Q. Li for shar- 
ing his unpublished data. This work was partially supported 
by the Strategic Japanese-Finland Cooperative Program on 
"Functional Materials," JST, Japan. 



1 G. Appli and Z. Fisk, Comments Condens. Mater. Phys. 16, 15 5 34 C. Herring, Phys. Rev. 96, 1163 (1954). 
(1992). 

2 M. Kasaya, F. Iga, M. Takigawa, and T. Kasuya, J. Magn. Magn. 
Mater. 47&48, 429 (1985). 

3 M. F. Hundley, P. C. Canfield, J. D. Thompson, Z. Fisk, and J. M. 
Lawrence, Phys. Rev. B 42, 6842 (1990). 

4 T. Takabatake, Y. Nakazawa, and M. Ishikawa, Jpn. J. Appl. Phys. 
26, Suppl. 547 (1987). 

5 T. Takabatake, F. Teshima, H. Fujii, S. Nishigori, T. Suzuki, T. 
Fujita, Y. Yamaguchi, J. Sakurai, and D. Jaccard, Phys. Rev. B 41, 
9607 (1990). 

6 V. Jaccarina, G. K. Wertheim, J. H. Wernick, L. R. Walker, and S. 
Arajs, Phys. Rev. 160, 476 (1967). 

7 D. Mandrus, J. L. Sarrao, A. Migliori, J. D. Thompson, and Z. 
Fisk, Phys. Rev. B 51, 4763 (1995). 

8 D. Zur, D. Menzel, I. Jursic, J. Schoenes, L. Patthey, M. Neef, K. 
Doll, and G. Zwicknagl, Phys. Rev. B 75, 165103 (2007). 

9 M. Klein, D. Zur, D. Menzel, J. Schoenes, K. Doll, J. Roder, and 
F. Reinert, Phys. Rev. Lett. 101, 046406 (2008). 

10 C. Petrovic, Y. Lee, T. Vogt, N. Dj. Lazarov, S. L. Bud'ko, and P. 
C. Canfield, Phys. Rev. B 72, 045103 (2005). 

11 A. Bentien, S. Johnsen, G. K. H. Madsen, B. B. Iversen, and F. 
Steglich, Europhys. Lett. 80, 17008 (2007). 

12 P. Sun, N. Oeschler, S. Johnsen, B. B. Iversen, and F. Steglich, 
Appl. Phys. Express 2, 091 102 (2009). 

13 P. Sun, N. Oeschler, S. Johnsen, B. B. Iversen, and F. Steglich, 
Phys. Rev. B 79, 153308 (2009). 

14 P. Sun, N. Oeschler, S. Johnsen, B. B. Iversen, and F. Steglich, J. 
Chem. Soc, Dalton Trans. 39, 1012 (2010). 

15 A. Perucchi, L. Degiorgi, R. Hu, C. Petrovic, and V. F. Mitrovic, 
Eur. Phys. J. B 54, 175 (2006). 

16 A. Herzog, M. Marutzky, J. Sichelschmidt, F. Steglich, S. Kimura, 
S. Johnsen, and B. B. Iversen, Phys. Rev. B 82, 245205 (2010). 

17 A. Bentien, G. K. H. Madsen, S. Johnsen, and B. B. Iversen, Phys. 
Rev. B 74,205105 (2006). 

18 H. Takahashi, Y. Yasui, I. Terasaki, and M. Sato, J. Phys. Soc. Jpn. 
80,054708 (2011). 

19 J. F. Woods and C. Y. Chen, Phys. Rev. 135, A1462 (1964). 

20 A. Kawabata, J. Phys. Soc. Jpn. 49, 628 (1980). 

21 N. W. Ashcroft and N. D. Mermin, Solid State Physics (Holt, 
Rinehart and Winston, New York, 1976). 

22 P. J. Moria and J. P. Maita, Phys. Rev. 96, 28 (1954). 

23 P. P. Debye and E. M. Conwell, Phys. Rev. 93, 693 (1954). 

24 T. Takabatake, Y. Echizen, T. Yoshino, K. Kobayashi, G. 
Nakamoto, H. Fujii, and M. Sera, Phys. Rev. B 59, 13878 (1999). 

25 M. Sasaki, G. X. Tai, S. Tamura, and M. Inoue, Phys. Rev. B 47, 
6216 (1993). 

26 R. Hu, V. F. Mitrovic, and C. Petrovic, Appl. Phys. Lett. 92, 
182108 (2008). 

27 J. S. Kim, J. Appl. Phys. 86, 3187 (1999). 

28 G. S. Nolas, J. Sharp, and H. J. Goldsmid, Thermoelectrics: Ba- 
sic Principles and New Materials Developments (Springer, Berlin, 
2001). " 

29 K. Wang, R. Hu, and C. Petrovic, larXiv: 1 102.05531 (to be pub- 
lished). 

30 J. M. Tomczak, K. Haule, T. Miyake, A. Georges, and G. Kotliar, 
Phys. Rev. B 82, 085 104 (2010). 

31 H. P. R. Frederikse and E. V. Mielczarek, Phys. Rev. 99, 1889 
(1955). 

32 H. P. R. Frederikse, Phys. Rev. 92, 248 (1953). 

33 L. Weber and E. Gmelin, Appl. Phys. A 53, 136 (1991). 



